Creation of Dirac monopoles in spinor Bose-Einstein condensates 
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We demonstrate theoretically that using standard external magnetic fields, one can imprint point-like topo- 
logical defects to the spin texture of a dilute Bose-Einstein condensate. Symmetries of the condensate order 
parameter render this topological defect to be accompanied with a vortex filament corresponding to the Dirac 
string of a magnetic monopole. The vorticity in the condensate coincides with the magnetic field of a magnetic 
monopole, providing an ideal analogue to the monopole studied by Dirac. 
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The existence of particles with non-zero magnetic charge, 
that is, magnetic monopoles has far-reaching implications to 
the laws of quantum mechanics, theories of elementary par- 
ticles, and cosmology [ll Q S- However, experimental ev- 
idence of magnetic monopoles as fundamental constituents 
of matter is still absent, and hence there is great incentive to 
search corresponding configurations in experimentally more 
tractable systems. Several aspects of monopoles have been 
investigated in the context of liquid crystals anomalous 
quantum-Hall effect [5], and exotic spin systems |6]. Re- 
cently, Qi et al. proposed that a magnetic monopole could 
be induced in a topological insulator |7]. Despite these el- 
egant proposals, an experimental realization of a magnetic 
monopole as an emergent particle or any analog of the Dirac 
monopole in still lacking. One of the most promising can- 
didate systems to realize such analogies has been superfluid 
•^He BSllslllOllllI], but to date there are no direct experimental 
observations of these topological excitations. 

Dilute Bose-Einstein condensates (BECs) of alkali atoms 
with a hyperfine spin degree of freedom combine magnetic 
and superfluid order and share many common features with 
the superfluid ^He | IJ, Y2, [Tsl [3. The order parameter 
describing such systems is typically invariant under global 
symmetries that form a non-Abelian group and monopoles 
and other textures can occur if this symmetry is broken. On 
the other hand, spinor BECs are well-suited to host artifi- 
cially generated gauge fields which can provide an alternative 
method to realize a magnetic monopole |15]. In the simplest 
case of spin-1 condensate, a variety of different topological 
defects such as global monopoles |16ll . non-Abelian magnetic 
monopoles [17], global textures 11811 and in particular, analo- 
gies to the Dirac monopole [ 191 have been investigated. In the 
related two-component condensates, skyrmion textures have 
been studied by several authors 1 2^ 21 1 22 1 . An experimen- 
tal realization of any of these topological states still remains a 
milestone in the field of cold atoms. 

In this Letter, we consider a spinor BEC with the total hy- 
perfine spin F = \ which in the absence of external magnetic 
fields has two phasesj^ a ferromagnetic and an antifeiTomag- 
netic (polar) phase il3ll . However, in the presence of a strong 
enough external magnetic field, the spin of the order parame- 



ter aligns with the local field and the condensate order param- 
eter corresponds to the ferromagnetic phase. Modifying the 
external field adiabatically, multi-quantum vortices can be im- 
printed into the condensate as was proposed theoretically by 
Nakahara et al. lEsll and verified experimentally by Leanhardt 
et al. ||24|[ . The method we employ here utilizes the same ide- 
ology in this respect, but due to non-trivial three-dimensional 
structure of the magnetic field we are able to create a point- 
like defect to the spin texture of the condensate that gives rise 
to vorticity coinciding with the magnetic field of a magnetic 
monopole. 

Let us assume first that the hyperfine spin of the condensate 
is aligned with an external magnetic field which is a combina- 
tion of two quadrupole fields and a homogeneous bias field 



B{r,t)=B[{xx + yy) 



B'^zz 



Bo{t)b, 
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where Maxwell's equations impose the condition 2B[ + B2 — 
for the magnetic field gradients. The direction of the 
bias field is determined by the unit vector b. Such combi- 
nation of quadrupole fields is produced by, e.g., a crossing 
pair of Helmholtz coils, and it has recently been proposed to 
dynamically generate knot-like textures in antifeiTomagnetic 
BECs llStl . The point where the external magnetic field van- 
ishes is the center of a monopole defect in the spin texture 
S — ^^^'^ given by the condensate order parameter * = 
(-01, tpo, tp-i)'^ and the spin-1 matrices ^ — {Tx, ^^j, ^z)- 
This monopole can be characterized by the charge 119(1 



(2) 



where s = S/\S\ and the integral is taken over a surface 
S enclosing the defect. The indices take values 1,2,3 and 
summation over repeated indices is implied. The Levi-Civita 
tensor is denoted by e^p^. For the magnetic field in Eq. ([T]), 
the charge of the spin texture becomes Q = 1. A schematic 
illustration of two possible textures is shown in Fig. [T] For 
the feiTomagnetic order parameter manifold, the second ho- 
motopy group is trivial and isolated monopole defects are not 
allowed. Thus the monopole defect in the spin texture is as- 
sociated with a vortex filament extending outwards from the 
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FIG. 1: Possible configurations for the spin texture. Both configura- 
tions are characterized by the same charge Q = 1 but only the one in 
(a) can be imprinted using the combination of two quadrupole fields. 
The hedgehog texture in (b) describes the vorticity fts correspond- 
ing to the spin texture in panel (a). In both cases, the vector fields 
are symmetric with respect to rotations about the axis on which the 
singularity filament (zigzag line) lies. 



monopole giving rise to a physical Dirac string fEllsl]. Since 
the monopole is located at the zero of the magnetic field, ad- 
justing the bias field Bq (t) moves the monopole along the axis 
given by b. In particular, the monopole can be brought in from 
outside of the condensate by taking large enough bias to begin 
with and then ramping it down adiabatically. 

The superfluid velocity can be defined as = 
—ih(^V(/m, where \I> = ^/qC, Q is the derisity of particles, 
and m is the mass of the constituent atoms lll3lll9l l251. For 
simplicity, we first assume that b ^ z and define new coordi- 
nates by {x',y',z') — {x,y,2z — Bo{t)/ B[). In the adiabatic 
limit, the condensate order parameter corresponds to the local 
eigenstate of the linear Zeeman operator gj^/i^ |B| n ■ T with 
h = (sin (3 cos a, sin (3 sin a, cos /3). The local Zeeman state 
can be written as 



|C(a,/3)) =Uia,p)\F,mt 



(3) 



where the unitary transformation U{a,/3) is given by 
U{a,l3) = e~'"-^'-e^*'3^''e*"^- and the state \F,mp) is an 
eigenstate of JF^ with an eigenvalue rup — —F, • • • , +F. 
Assuming that the local Zeeman energy is minimized for 
B[ < Q, Bq > 0, and gp < Owe obtain in the F — 1 case 



n 1 -cosi9' 
m r' sin z9' 



(4) 



where {r',(p','d') refer to spherical coordinates in the new 
coordinate system. Similar result has also been obtained in 
Ref. | | 19.1 for the hedgehog texture of Fig. [Tib). 

The superfluid velocity is equivalent to the vector poten- 
tial of a magnetic monopole and it has a Dirac string along the 
negative z' axis if Bo{t = 0) > 0. For spinor condensates the 
superfluid flow can be irrotational and it is characterized by its 
vorticity fJ^, = V x which becomes 



indicating that vorticity corresponding to the imprinted 
monopole defect is equivalent to the magnetic field of a mag- 
netic monopole, see Fig.[T] In particular, the topology of Jl^ 
is unaffected by the scaling and translation, and hence Jl^ re- 
mains equivalent to the hedgehog texture shown in Fig. [Tib) 
also in the original coordinate system. The Mermin-Ho rela- 
tion for general spin-F 1 25 , 2^ can be computed from Eq. ([3]) 
yielding 
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In the adiabatic limit spin S aligns with h and the charge Q in 
Eq. (O is thus directly proportional to the flux of il^ through 
a surface E enclosing the defect. The flux of the monopole is 
supplied by the Dirac string which is omitted form Eq. ((5]l. In 
case of F = 1, the monopole flux is 2h/m, that is, two angu- 
lar momentum quanta, which implies that the vortex filament 
terminating at the monopole must carry the same amount of 
vorticity. The situation is thus similar to the Dirac monopole 
in superfluid '^He-A For a general spin-F BEC the 

Dirac string carries 2F quanta of angular momentum. 

Non-adiabatic effects arising from interactions between 
atoms, kinetic energy, and the finite timescales in manipulat- 
ing the external magnetic fields can render the spin to deviate 
from the direction of the local magnetic field. We take these 
effects into account by solving the dynamics of the spinor 
order parameter from the Gross-Pitaevskii mean-field equa- 
tion 01 [il 



0„ 



h 1 



(5) 



ift^*- [/io + Co|*|2 + c2(*t:P*).:P]*, (7) 

where the interaction strengths Cq and C2 depend on the scat- 
tering lengths in the different channels coiTesponding to the 
total hyperfine spin of two scattering particles 1. 1 3,1 . The 
single-particle operator is given by 

ho = -T^V^ + Vopt + g^-fisB ■ 
2m 

and in the simulation, we employ an external optical poten- 
tial Vopt = mr'^Ld'^ /2 and choose the parameters according to 
*^Rb which implies ferromagnetic interactions. 

In the simulation, energy time, and spatial variables are 
given in the units of hw^., l/w^ and = ^Jh/muj^., re- 
spectively. For = 27r X 250, the dimensional values of 
the parameters are given by B'^ — —0.05 T/m and -Bo(i = 
0) = 1 |jT corresponding to 10'' atoms. We have also con- 
sidered a different atom number and antiferromagnetic in- 
teractions [l33ll . The area considered in the simulation is 
24 X 24 X 27 in the units of and the size of the compu- 
tational grid varies from 141 x 141 x 161 to 175 x 175 x 195 
points. In the simulation, we first calculate the ground state of 
the system corresponding to the initial values of the magnetic 
fields using the successive over relaxation (SOR) algorithm, 
and then propagate the initial state according to the time- 
dependent GP equation ((7]i using the split operator method 
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FIG. 2: Particle density at different stages of monopole creation. The 
monopole is created by ramping tlie bias field down in a period to — 
SO/tJr. Blue (dark) color denotes densities from g = 3.8 x 10~^ 
N/a^ to g = 1.8 X 10"'' N/a^ and green (light) color densities 
from g ^ lAx 10"^ N/a^ to g = 1.5 x 10"^ N/af.. The Dirac 
string is identified as the density depletion that propagates through 
the condensate. 



combined with the Crank-Nicolson method. The time step 
used in the simulation is 10~'^/uj^. 

Since the Dirac string carries two quanta of angular mo- 
mentum it is expected to be prone to splitting into two separate 
strings each carrying one angular momentum quantum. To 
avoid this scenario, we first imprint the monopole defect with 
the bias field parallel to z axis. The particle density from this 
simulation is shown in Fig.|2] The bias field is ramped linearly 
down to zero in a period tQ = 50 /ujr and then equally to neg- 
ative values. The spin texture <S takes the form described in 
Fig- Ufa) and the corresponding vorticity is shown in Fig.|3la). 
Since the monopole is brought in along the symmetry axis 
of the system, the Dirac string remains intact. The created 
monopole appears to be relatively stable: after the bias field 
was ramped down, the monopole was allowed to evolve in 
time for a period in the presence of the two quadrupole 
fields resulting in slow small-amplitude oscillations along the 
z axis but it remained otherwise intact. Due to nonadiabatic 
effects in the creation process, the monopole defect in S lags 
slightly behind the zero point of the magnetic field and the 
Dirac string deviates from a pure 5-function distribution. We 
have also taken the parameters according to ^^Na, i.e., corre- 
sponding to antiferromagnetic interactions, and obtained es- 
sentially the same behavior as in the case of ^^Rb. 

If the rotation symmetry with respect to z axis is explicitly 
broken by a bias field with b x £ 7^ 0, we observe the Dirac 
string to split into two parts, see Fig|3jb). The particle den- 
sity is not depleted all the way along the two strings suggest- 
ing that nonadiabatic effects become more pronounced in this 
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FIG. 3: Vorticity corresponding to the monopole. Vorticity fl^ is 
computed numerically after the bias field is ramped from Bq (0) = 
6 . 7 /iT to _Bo (io ) ~ and it is shown in the units oih/ma^. (a) to = 
50/LJr and the bias field is parallel to z axis, (b) to — 40 /ui^ and the 
bias field is parallel to the vector 3/2 i + y + 14/3 £. Magnitude 
of the vorticity is denoted by color and the map is linear between the 
minimum and maximum values. For clarity, only the relevant parts 
of the vector field f2, are shown. 



case. The spin density on the contrary, exhibits a clear 
depletion along the path of the two Dirac strings [32J and to- 
gether with the state-of-the-art experimental methods 1 3^ 31 ] 
it should give an efficient signature of the monopole. In the 
experiments, a typical imperfection is a slight misalignment 
and slow drift between the center of the optical trap and the 
symmetry axis of the magnetic field. We model this by tak- 
ing b\\z and adding a small constant term Bq x to the mag- 
netic field in Eq. ([T]!. This introduces an offset equal to 7% 
of the effective radius of the condensate to the path traced by 
the monopole with respect to z axis. The offset breaks the 
rotation symmetry of the system but the Dirac string of this 
off-axis monopole remains undivided (data not shown) sug- 
gesting that the instability of the string becomes visible only 
for fairly large perturbations. 

Since the monopole defect in the spin texture is topologi- 
cally unstable, it can be removed by local surgery. Thus we 
address to the fate of the monopole after it has been created 
and all the magnetic fields pinning the monopole are turned 
off. In the case of superfluid ■^He, Dirac monopoles are ex- 
pected to drift to the boundary of the vessel to form a surface 
defect — ^boojum |27]. To analyze the decay of the monopole 
in our case, we have carried out simulations in which the ex- 
ternal magnetic fields are switched off immediately or ramped 
down with constant speed in a period <i. For the finite switch- 
off times, the decay of the monopole initiates already while 
the external fields are being ramped down. The initial state in 
the simulation corresponds to a monopole created by ramping 
the bias field aligned with z axis down in a period to = 50/w^. 

For ferromagnetic interactions, we have used three differ- 
ent values ti = 0,ti — 2.5/ojj. and ti = 5.0/a;^. In all three 
cases, the qualitative features are the same and they are shown 
schematically in Fig.|4](see also 13211 ). For the monopole de- 
fect in the spin texture, the monopole unwinds itself along the 
Dirac string and results in a closed vortex ring, see Fig.|4|a-d). 



4 



In the course of unwinding, a cylindrical domain wall separat- 
ing the expanding core of the Dirac string from the rest of the 
texture is formed. Depending on how fast the external fields 
are turned off, the domain wall is either directly pushed out 
of the condensate (slow turn-off) or contracted to another vor- 
tex ring which eventually drifts out of the condensate (rapid 
turn-off). In the simulations, the resulting vortex ring in spin 
texture persists until the end of the simulation which in all 
three cases spans a period lO/o;^. On the other hand, the un- 
winding of the monopole in the vorticity is depicted in 
Fig.Ue-h). During the unwinding of the monopole, vorticity 
concentrated at the Dirac string diffuses outwards and tends 
to relax towards more uniformly distributed values. 



Einstein condensates. Our proposal is directly realizable 
in the present-day experiments and recent developments in 
the non-destructive imaging of the magnetization of spinor 
BECs isoi 31 1 pave way for detecting the exciting dynamics 
of these monopoles such as the unwinding of the monopole. 
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FIG. 4: Unwinding of the monopole defect in the spin texture S (a)- 
(d) and in the vorticity fig (e)-(h). Arrows pointing upwards (down- 
wards) from the plane are denoted by (0). Along the dashed line 
in (b) and (c), the spin points either upwards or downwards from 
the plane. The unwinding in fig is shown in a plane perpendicu- 
lar to 2 axis and it is qualitatively independent of the z coordinate. 
High-resolution figures from the simulation corresponding to these 
schematic illustrations are available in [32]. 

For antiferromagnetic interactions we have carried out a 
simulation with magnetic fields ramped down linearly in a pe- 
riod ti = 2.5 /u!^. The results agree qualitatively with the 
ferromagnetic case and the end configuration is again a vor- 
tex ring in the spin texture. Another vortex ring is gener- 
ated at the boundary of the condensate and it persists until 
the end of the simulation. The length of the simulation was 
in this case. It should be noted that simulations con- 
cerning the unwinding of the monopole span a relatively short 
period and thus it remains a question for the future research to 
study if the evolution of the created monopoles differs qual- 
itatively between feiTomagnetic and antifeiTomagnetic inter- 
actions. From the global symmetries of the order parameter 
in the absence of external fields ITs, 28 1, one could expect 
different behavior by purely topological reasoning |4]. Fur- 
thermore, we have not included the possible magnetic dipole- 
dipole interaction which may change the dynamics of the 
monopole unwinding. In particular, the spin texture consid- 
ered here resembles the co-called two-z-flare texture that was 
found to be stabilized by the dipole-dipole interaction |29]. 

In conclusion, we have introduced and modelled a ro- 
bust method to create Dirac monopoles in spinor Bose- 
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t = 38/ w,, t = 40 /u!,. 

FIG. 1: Spin density of the condensate. Splitting of the Dirac string is manifested in the spin density \S\ which is depleted along the two 
strings. The figure corresponds to spin densities from \S\ — 1.5 x 10^* N/a^ to |5| = 4.0 x 10^* N/a^. Here x axis points downwards 
from the plane. 
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FIG. 2: Unwinding of tlie monopole defect in tlie spin texture. The 2D plane corresponds to x = 0. In the figure, Sy and Sz components of 
the unit spin vector s = are shown. The spatial scale in different panels changes since the radius of the vortex ring changes during its 

formation. Time corresponding to each panel is 50.5/0;^ (a), (b), SS/cj^ (c), and GO/w^ (d). 
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FIG. 3: Unwinding of the Dirac monopole in vorticity. The 2D plane corresponds to ^ = — 2.8 a^. Vorticity is given in the units of h/ma^. 
The time instant corresponding to each panel is SO.S/oj^ (a), 51.5/0;^ (b), 52.5/cj,, (c), and 53.5/0;^ (d). At the center of the defect, the vector 
field is always parallel to the z axis pointing upwards from the plane. In the figure x and y axes correspond to horizontal and vertical axes, 
respectively. The field of view is ISa^ x ISa^. 



